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ABSTRACT 

Lynds dark cloud LDN1622 represents one of the best examples of anomalous dust 
emission, possibly originating from small spinning dust grains. We present Cosmic 
Background Imager (CBI) 31 GHz data of LDN1621, a diffuse dark cloud to the north 
of LDN1622 in a region known as Orion East. A broken ring with diameter 20 arcmin 
of diffuse emission is detected at 31 GHz, at « 20 — 30 mJy beam~^ with an angular 
resolution of w 5 arcmin. The ring-like structure is highly correlated with Far Infra- 
Red emission at 12 — 100 /im with correlation coefficients of r « 0.7 — 0.8, significant 
at ^ lOcr. Multi-frequency data are used to place constraints on other components of 
emission that could be contributing to the 31 GHz flux. An analysis of the GB6 survey 
maps at 4.85 GHz yields a 3cr upper limit on free-free emission of 7.2 mJy beam^^ 
(< 30 per cent of the observed flux) at the CBI resolution. The bulk of the 31 GHz 
flux therefore appears to be mostly due to dust radiation. Aperture photometry, at an 
angular resolution of 13 arcmin and with an aperture of diameter 30 arcmin, allowed 
the use of IRAS maps and the WMAP 5-year W-band map at 93.5 GHz. A single 
modified blackbody model was fitted to the data to estimate the contribution from 
thermal dust, which amounts to ~ 10 per cent at 31 GHz. In this model, an excess 
of 1.52 ± 0.66 Jy (2.3cr) is seen at 31 GHz. Future high frequency ~ 100 - 1000 GHz 
data, such as those from the Planck satellite, are required to accurately determine 
the thermal dust contribution at 31 GHz. Correlations with the IRAS 100 /im gave a 
coupling coefficient of 18.1 ±4.4 /iK (MJy/sr)~^, consistent with the values found for 
LDN1622. 

Key words: radio continuum: ISM - diffuse radiation - radiation mechanisms: gen- 
eral 



1 INTRODUCTION 

Anomalous Microwave Emission (AME) is the name given 
to excess microwave emission, observed at frequencies in 
the range ~ 10 — 60 GHz, that is strongly correlated 



with far infrared ( FIR) dust emissio n ( Leitch et al. 19971: 



Kogut et al.' 'l99^; 'Finkbcincr ct al' 2002; Banday et al 
'2003.: de Oliveira-Costa et al. 2004; Finkbeincr et al,, .20041:^ 
Watson et al. 2005: Casassus et al.l 



2008 



^^Davies et al.l 



20061: iDic kinson e t al.|]2009al : I Ami Consortium et aLlbOOa 
Tibbs et al. 2010). This dust-correlated emission is known 
to be a significant source of contamination for Cosmic 



Microwave Background (CMB) data, that must be sep- 
arate d accurately from t h e CMB signal IIBennett et aU 
20031: iBonaldi et al.1 [2OO7I: iMiville-Deschenes et all l20o'8r 
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Gold et al.ll2009l : IDickinson et all l2009bl ) . Ahhough there 
is still some debate about the physical mechanism that 
is responsible for the emission, the most favoured ex- 
planati on is in terms of small , rapidly spinning dust 
grains (jPraine fc LazarianI Il998al fbl) . Assuming this is the 
case, microwave observations of spinning dust represent a 
new way of studying the properties of interstellar dust 
grains and its environment within the interst e llar medium 
llAli-Haimoud et al.1 l2009l : iDobler et all l2009l : lYsard et all 
l201Cf ). Accurate observations, at frequencies in the range 
5 — 100 GHz, are now required to confirm spinning dust 
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grains as the source of the anomalous emission and to study 
the spectrum to infer information about the grains and their 
environs. 

One of the best examples of spinning dust emis- 
sion, comes from observations of the dark c loud LDN1622 
ijFinkbeiner et al.ll2002l : ICasassus et al.ll2006h . LDN1622 is a 
small (« 10 arcmin) dark cloud at the low Galactic latitude 
edge of the giant m olecular cloud Orion B, o n the outer edge 
of Barnard's loop (|Maddalena et aDll986l ). The spectrum 
between 1 GHz and 3000 GHz is well-fitted by a superpo- 
sition of optically thin free-free emission, thermal dust at a 
temperature of T ~ 15 K and a spinning dus t component 
with a peak at ~ 30 GHz (Casa ssus et al.|[2006l ): at 30 GHz, 
and on angular scales < 20 arcmin, the spectrum is domi- 
nated by spinning dust emission. Furthermore, correlations 
with IRAS infra-red maps indicate a better correlation with 
the emission from Very Small Grains (VSGs), as expected if 
the origin of the cm emission is spinning dust. 

About 25 arcmin to the north-east of LDN1622 lies 
another dark cloud, LDN1621, which is more diffuse 
than LDN1622 and forms a broken ring-like structure. 
These two clouds together were termed as Orion East by 
iHerbig fc Kameswara Raol (|l972l ) and are thought to lie 
at a distance of ~ 400 pc, although there is some de- 
bate about whethe r it lies much clo ser to us at ~ 140 pc 
jWilson et al.ll2005l : iKun et al.ll2008l '). LDN1622 is a higher 
density dust cloud in which low mass star formation is 
beginning, as indicated by the presence of T-Tauri stars 
iHerbig fc Kameswara R"aol 1 19721 ) but there are no bright 
OB stars in the vicinity; the UV radiation is predominately 
from the Orion OBI association. By contrast, LDN1621 is 
a lower density, ring-l ike cloud, with no obvious pre-main 
sequence stars within l|Lee et al.ll2005l ). 

Here we present 31 GHz data of the area around 
LDN1621 and compare it with multi-frequency data, to esti- 
mate the contributions of free- free, thermal dust and anoma- 
lous dust components. This is compared with previous obser- 
vations of the LDN1622 dark cloud to the south-west. Sec- 
tion [2] describes the observations and data reduction while 
Section [3] discusses multi-frequency maps of the LDN1621 
region. Section [4] gives a quantitative analysis of the contri- 
butions from free-free and thermal dust emissions and the 
possible origin of the 31 GHz excess. Conclusions are given 
in Section [5l 



2 OBSERVATIONS AND DATA REDUCTION 



The Cosmic Background Imager (CBI: iPadin et all (|2002l ')') 
is a 13-element close-packed interferometer, that operated at 
the Chajnantor Observatory, in the Atacama desert (Chile), 
from 1999-2008. The CBI was designed to image the CMB 
anisotropics on scales ~ 5 — 30 arcmin, in ten 1 GHz bands at 
26 — 36 GHz. With its high surface brightness sensitivity and 
calibration stability, it has also been used to image Galactic 
objects including diffuse dust clouds, supernova rernnants , 
planetary nebulae an d HII regions ^Casassus et al.' '2004 
2006 , 2007, 2008; Hal es et all |2004| ; [Dickinson et al . 2006; 
20071 . l2009al ). In 2006-2007, the CBI was fitted with 1.4 m 
dishes (CBI2), to allow optimal sensitivity on the smallest 
angular scales accessible by CBI (Taylor et al., in prep.). 
The larger dishes improved the point-source sensitivity and 



surface brightness sensitivity over angular scales ~ 5 — 20 ar- 
cmin. The primary beam is well-modelled by a Gaussian 
with a FWHM= 28.2 arcmin at 31 GHz up to ~FWHM/2, 
the region used in the analysis. 

During 2007 and 2008, CBI2 observed LDN1621 (J2000: 
R.A. 05''55™21.6' Dec. -h02°ll™33') over 6 nights with a 
combined integration time of 5.7 hours on source. A trail 
field, separated by 8 min in R.A., was observed at the same 
hour angles to allow removal of local correlated signals such 
as ground spillover. The data analysis follows a similar pro- 
ccdurc to that used in CMB data analysis fsee lSievers et all 
(2009) and references therein). The amplitudes and phases 
were calibrated using short observations of Tau-A, assum- 
ing a flux density of 341 Jy at 31 GHz and a spectral index 
a — —0.299 {S oc v"). The absolute calibration is known to 
0.5 per cent, with these values tied to the brightness tem- 
erature of Jup iter, assumed to be Tj = 146.6 ± 0.75 K 
Hill et al.|[2009l '). AU the data were inspected by eye to re- 
move noisy data and baselines/antennas that were not pro- 
viding adequate calibration. We assign an overall calibration 
uncertainty of 5 per cent to account for residual gain varia- 
tions, ground spillover and pointing errors. 

The calibrated visibilities were inverted to for m an im- 
age, a nd deconvolved using the CLE AN algorithm ( HogbomI 
Il974h . using the difmap program (|Shepherdlll997l V Fig. [T] 
shows the final CLEANed 31 GHz map of the LDN1621 
region. We chose natural weighting of the visibilities to op- 
timize for sensitivity rather than angular resolution. The 
CLEAN components were primary-beam corrected assum- 
ing a Gaussian model with FWHM= 28.2(!/ghz/31) ar- 
cmin. A low gain of 0.03 was used to recover as much 
extended emission as possible. The synthesized beam is 
6.1 X 4.8 arcmin and the noise r.m.s. level in the map is 
~ 6 mjy beam"^. Spectral indices within the CBI band 
were not reliable due to the extended nature of LDN1621 
and limited signal-to-noise ratio with a small frequency lever 
arm of 10 GHz. 



3 DISCUSSION AND COMPARISON WITH 
MULTI-FREQUENCY DATA 

The 31 GHz CBI image (Fig. [T]) shows a ring of emission 
associated with the dark cloud LDN1621. The ring, of di- 
ameter « 25 arcmin, is not complete and is not of uni- 
form brightness. The ring has a brightness in the range 
« 20 — 30 mJy beam~^, except for the northern part which 
is at ~ 35 mJy beam~^. It should be noted that the ring is 
not completely within the FWHM of the primary beam at 
31 GHz, as indicated in Fig.[T] The northern part of the ring 
is at « 15 arcmin from the pointing centre and is therefore 
likely to be more affected by primary beam errors and noise. 
Indeed, the brightest pixel in the image is 16 arcmin from 
the pointing centre. We found that too many iterations of 
the CLEAN algorithm could result in additional flux in pix- 
els outside the primary beam FWHM. In particular, for the 
brightest pixel, we found that the peak brightness could be 
boosted by up to a factor of ~ 1.7 compared to that pre- 
sented here. Also, far outside the primary beam, very bright 
sources could be artificially produced. For this reason, bright 
pixels that were significantly outside of the primary beam 
FWHM were not corrected for the primary beam, to prevent 
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Figure 1. CLEANed 31 GHz map of the LDN1621 region. Con- 
tours are at 40, 60 and 80 per cent of the peak brightness 
(35.4 mjy beam~^). The synthesized beam (6.1 X 4.8 arcmin) is 
shown in the lower-left corner. The FWHM of the primary beam 
is indicated by the thick black dashed line. Thin white solid lines 
indicate the apertures adopted in Section |4l 

distortion of the image. Inside the FWHM of the primary 
beam, such errors are neghgible. We also verified the mor- 
phology of the 31 GHz image with an image created using 
earlier CBIl data, of just 30 mins on a single night. Al- 
though of limited signal-to-noise ratio, the ring of LDN1621 
was detected with a similar morphology to the data pre- 
sented here, and the entire ring fits well within the primary 
beam FWHM of 45.2 arcmin at 31 GHz. It is also worth 
noting that LDN1621 is in fact visible on the edge of the 
CBI 3 1 GHz image of LDN1622 presented in' Casassus et all 
lj2006l ) and is seen to have a similar morphology. 

Fig. [2] shows a comparison of the 31 GHz contours 
(from Fig. [T| with various multi-frequency maps. The fre- 
quencies, angular resolutions and references are given in Ta- 
ble [1] The best IR data for comparison are the IRIS maps, 
which ar e the re-processed IRAS data at 12, 25, 60 and 
100 fim (jMiville-Deschenes fc Lagachell2005l ). We searched 
the Spitzer data archive for higher resolution IR data, but 
IRAC/MIPS data of the LDN1622 region only covered up 
to Dec.« -1-02° 10"" and therefore did not cover LDN1621. 
A ring of emission, coincident with the ring observed at 
31 GHz, is clearly seen in the IRIS maps at 12—100 /im. The 
IR ring of LDN1621 has a similar morphology to the 31 GHz 
contours. In particular, the ring is not complete to the west 
and has enhanced brightness on the north and south sides. 
The bright cloud to the south-west of the ring is the dark 
cloud LDN1622, and is not detected in these CBI data be- 
cause of the primary beam attenuation beyond « 20 arcmin 
from the pointing centre. 

The low frequency radio data at 1.4, 2.3 and 5 GHz 
(Fig. [21 do not show any evidence of emission associ- 
ated with LDN1621. The relatively high angular resolution 
(45 arcsec) of the 1.4 GHz data is ideal for identifying ex- 



tragalactic radio sources in the field that could be confusing 
the extended Galactic emission. No bright radio sources are 
detected in the field and no radio sources are coincident with 
the hot spots see n at 31 GHz down t o a level of 2.5 mJy in 
the NVSS maps (ICondon et aLlliggsl ). The majority of ex- 
tragalactic sources have steeply falling spectra and therefore 
are unlikely to be detected at 31 GHz. Due to the limited 
u, V coverage of the NVSS data, extended emission on scales 
of a few arcmin and larger will be resolved out by the in- 
terferometric response of the VLA, particularly in snap-shot 
mode. However, the 4.85 GHz GB6 image, at a resolution 
of 3.5 arcmin, also shows no evidence of emission in the 
vicinity at the level of > 8 mJy beam"'^. The GB6 maps 
are also affected by spatial filtering on scales > 20 arcmin 
l|Condon et al.l [19941 ) but remain sensitive to spatial scales 
comparable to the CBI. 

The WMAP 5-yr maps are the only other data available 
at frequencies close to 31 GHz. WMAP maps at Ka-band 
(33 GHz) and Q-band (40.7 GHz) show enhanced emission 
on large angular scales, while the ring is barely visible be- 
cause of the limited angular resolution of the WMAP data 
(Table [T}. At W-band (93.5 GHz), the large scale emission 
becomes faint and only LDN1622 is detected. In Section r4.2l 
we use the W-band data to constrain the contribution from 
thermal dust. 

The optical images of the LDN1621 region (Fig. ^ 
clearly show extended emission in the vicinity of LDN1621 
and LDN1622. Most of this will be emission from the 
bright Ha line, which has been separated in the SHASSA 
continuum-subtracted Ha map. The Ha line is a good tracer 
of warm ionized gas and because of its well-known depen- 
dence on the electron dens ity (Jhq oc fn^dl), it is a good 
tracer of free-free emission l|Dickinson et al.ll2003l '). There is 
a strong gradient of large-scale emission, becoming stronger 
to the south-west. This is expected since there is a strong 
UV-radiation field towards the Orion OBI association in this 
direction. A large-scale gradient such as this will not be de- 
tected by the CBI as it is not sensitive to scales > 20 arcmin. 
An enhancement of Ha intensity is visible around LDN1622 
while it appears to be absorbing the majority of Ha photons 
from LDN1622, which suggests that LDN1622 is in front of 
the ionized gas responsible for the Ha emission. A similar 
situation occurs in the ring of LDN1621 but the effect is 
much less pronounced because the Ha intensity is less and 
there is a stronger gradient at the position of LDN1621. 
Indeed, some estimates of the distance to LDN1621/1622 
put it as close as 1 20 pc compared to 400 pc for Orion B 
(| Wilson et al.ll2005^ . 



4 ANALYSIS & DISCUSSION 
4.1 Limits on free-free emission 

The GB6 4.85 GHz data provide a strong upper limit on the 
contribution of free-free emission at 31 GHz. Steeper spec- 
trum emission such as synchrotron radiation will be negligi- 
ble at 31 GHz. No emission is detected within the LDN1621 
area. After smoothing to the CBI resolution, the r.m.s. noise 
level at 4.85 GHz is « 3 mJy beam~^ corresponding to a 3a 
upper limit of 9 mJy beam"'^. At frequencies > 1 GHz, free- 
free emission will be optically thin, except for the densest 
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Figure 2. Multi-frequency maps of the LDN1621 region, superposed with the CBI 31 GHz contours from Fig. [T] The maps are, with 
colour scale units in parentheses, as follows: NVSS 1.4 GHz (mjy beam~^), HartRAO 2.3 GHz (mK), GB6 4.85 GHz (mjy beam~^), 
SHASSA continuum-subtracted H« (R), DSS red (arbitrary), WMAP Ka-band (/iK), WMAP Q-band [pK), WMAP W-band IRIS 
12 /im (MJy/sr), IRIS 25 /^m (MJy/sr), IRIS 60 ^^^a (MJy/sr) and IRIS 100 ^^ra (MJy/sr). Each map is at the original angular resolution 
(see Table [TJ. 



Table 1. Ancillary da ta used 
fl]: ICondon et al.l l|l998l l: 



in thi s paper. R ef erence s 
[2]: I Jonas et al.l lll99Sl) 



[3]: ICondon et al.l lll994l'): [4]: iHinshaw et al.l ll200g|): [5] 

iMivil lc-Dcschcncs fc Lag achj boOa) : [6]: iGaustad et al.l I2OO1I) 
[71: [McLe an et a l., (2000). 



Dataset 


Frequency/ 


Bcamwidth 


Reference 




Wavelength 


(FWHM arcmin) 




NVSS 


1.4 GHz 


0.75 


[1] 


HartRAO 


2.3 GHz 


20 


[2] 


GB6 


4.85 GHz 


3.5 


[3] 


WMAP Ka 


33.0 GHz 


« 37 


[4] 


WMAP Q 


40.7 GHz 


fs 29 


[4] 


WMAP W 


93.5 GHz 


« 13 


[4] 


IRIS12 


12 /^m 


3.8 


[5] 


IRIS25 


25 /^m 


3.8 


[5] 


IRIS60 


60 /xm 


4.0 


[5] 


IRISlOO 


100 


4.3 


[5] 


Ha 


656.2 nm 


Ri 4 


[6] 


DSS Red 


571 nm 


« 0.03 


[7] 



clouds, such as ultracompact HIT regions. We can therefore 
reliably extrapolate the upper limit to 31 GHz assuming a 
flux density spectral index {S oc of a = —0.12, appro- 
priat e for optically thin free-free emission (jPickinson et al.l 
l2003h . The 3ct upper limit to free-free emission at 31 GHz 
is then 7.2 mJy beam^^. We can therefore be quite confi- 
dent that the bulk of the 31 GHz emission cannot be due 
to synchrotron or free-free emission, although a small (up 
to a maximum of ~ 30 per cent) contribution from free-free 
contribution cannot be ruled out with the data available. 

The level of free-free emission can also be estimated 
from Ha data. If we assume that all of the absorbing dust 
lies in front of LDN1621, we can estimate the true Ha in- 
tensity by estimating the dust extinction along this line of 
sight. Along the LDN1621 ring, the total Ha intensity is 
/hq ~ 40 — 70 R. This is below the background value of 
~ 100 R and therefore the CBI will only detect a small 
fraction of this on scales « 4 — 20 arcmin. The dust absorp- 
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tion can b e estimated using the column-density map, at 
100 ^m, of ISchlegel etH] j 19981 ) assuming typical ISM con- 
ditions and reddening values. If we assume all the absorbing 
dust is in front of the Ha-emitting gas, the corrected Hq in- 
tensity is given by /g'^" = /hc x IqO-O^^'''^'^ (jPickinson et all 
I2OO3I ). where is the temperature-corrected intensity map 
in units of MJy/sr. Along the ring, D'^ ~ 50 - 60 MJy/sr, 
thus I^^'' ~ 500 R. For a typical electron temperature of 
7000 K, we expect a 31 GHz br ightness temperature of 
5.44 fiK/R (|Dickinson et al.ll2003l '). or 2.7 mK at 31 GHz. 
This corresponds to ~ 200 mjy beam"^ at 31 GHz at the 
CBI resolution. However, this estimate is based on assuming 
all the dust along the line-of-sight is absorbing and that all 
the Ha-correlated emission is detected by the CBI. We can 
see from the SHASSA image that the bulk of Ha emission is 
on large-scales and thus would not be detected by the CBI. 
The fluctuations in Ha on scales of ~ 5 arcmin are at the 
10 — 20 R level, which corresponds to ~ 10 — 20 mJy beam~^ 
at 31 GHz, if all the dust is in front of the Ha-emitting gas. 

We must therefore rely on the low frequency data for the 
strongest constraints on free-free emission on these angular 
scales. We also note that there is a filament of Ha that enters 
the the LDN1621 ring at a level of ~ 100 R (see Fig. ^ but 
is not detected by the CBI. Given the noise level in the 
31 GHz image is ~ 6 mJy beam~^, this could be considered 
an approximate upper limit on the the free-free contribution, 
which is clearly small on CBI angular scales. 



4.2 Limits on vibrational (thermal) dust emission 

The 31 GHz data are highly correlated with the IR data at 
12 — 100 fim. The IR data are dominated by thermal emis- 
sion from dust grains heated by UV radiation. The emission 
is often approximated by a modified blackbody function. 
Id oc h'^^'^^B{i/,Td), defined by an emissivity index, /3, and 
a blackbody dust temperature, T^. With typical ISM dust 
temperatures in the range Td ~ 10 — 50 K, there is a peak at 
~ 3000 GHz (100 Aim). Typically, the Ray leigh- Jeans tail of 
this emission dominates down to ~ 100 GHz and becomes 
negligible at lower frequencies where synchrotron and free- 
free emission dominate. For dust grains in the ISM, an emis- 
sivity index of /3 ~ 2 is typical, although a range of values has 
been observed, reflecting different grain properties and size 
populations (see , Pupae et al. ( 2003 ) and references therein). 
It is therefore conceivable that the CBI 31 GHz data could 
simply be the long wavelength tail of the thermal dust. 

With no high resolution IR data available below 
3000 GHz, it is somewhat difficult to extrapolate accurately 
the thermal dust from IRAS data alone; we would need to 
assume a fixed value for the emissivity inde x. However, the 
highe st frequency channel of WMAP data l|Hinshaw et al.l 
bOOgl ) at W-band (93.5 GHz) has a resolution of « 13 ar- 
cmin and will be useful for placing limits on the level of the 
thermal dust contribution, at least on these angular scales. 
Fig. [1 shows the WMAP 5-year 93.5 GHz map. There is 
no strong emission detected at the location of LDN1621, al- 
though LDN1622 is visible to the south-west and there is 
low level emission just to the south of LDN1621. 

To make a more quantitative estimate of the thermal 
dust contribution, we smoothed the IRIS maps to a com- 
mon resolution of 13 arcmin to allow a comparison with the 



93.5 GHz data. We then measured the flux density in an 
aperture of diameter 30 arcmin centred on LDN1621 (see 
Fig. H}. An estimate of the background level was estimated 
from an equivalent aperture placed in the top-right hand 
corner of the image (R.A. 05''53™14= Dec. -h02°33'"18') and 
subtracted from the LDN1621 aperture. Errors were derived 
from the r.m.s. fluctuations in the map. The results are given 
in Table [2l To obtain an equivalent 31 GHz value, we need 
to account for flux that has been lost due to the incom- 
plete u, u-coverage of the CBI observations. We therefore 
estimated the flux density based on the map at 100 fim, 
scaled to 31 GHz via the correlation coefficient derived in 
Section 1431 of 18.1 ±4.0 fiK (MJy/sr)"^; note that this is 
equivalent to calculating the amount of flux loss based on 
the 100 /im morphology and correcting the CBI observed 
integrated flux density. The aperture analysis yielded an in- 
tegrated flux density of 1.67 ±0.66 Jy at 31 GHz. The error 
includes the r.m.s. from the sky and the error in the corre- 
lation coefficient. 

The spectrum based on the aperture analysis is shown 
in Fig. O We fitted a modified blackbody curve to the 93.5, 
2997 and 4995 GHz data points only. The higher frequency 
data points (12/25 fim) do not give a good fit to such a sim- 
ple model because the grains are not in thermal equilibrium 
with the radiation field and there are emission lines (ionic 
and PAH) within the 12/25 /im passbands of IRAS. We do 
not consider the free-free component here because the GB6 
4.85 GHz data are not reliable when integrating over regions 
> 20 arcmin. The best-fitting model provides a good fit to 
the data with Td = 38.5 ± 4.3 K and /3 = 2.23 ± 0.18. The 
dust temperature is somewhat higher than expected for a 
dark cloud; we may be seeing the effects of the local diffuse 
ISM in the region which may be warmer than the dust grains 
truly associated with LDN1621. The coupling coefficient is 
close to the expected value, and is constrained mostly by the 
W-band data point. We also tried fixing the coupling coeffi- 
cient and/or emissivity index, to see if an adequate fit could 
be found but with a flatter index. However, it is not possi- 
ble to find a different best fit to all three data points with 
just a single dust component. An example is shown in Fig. [3] 
where the emissivity index was fixed at /3 — ±1.7; the best- 
fit is clearly not consistent with the data and the excess at 
31 GHz still remains. However, we do note that multiple dust 
components could, in principle, still contribute at 31 GHz. 
An example would be a colder (Td < 10 K) component, yet 
it cannot be dominant since the total emission must remain 
consistent with the 93.5 GHz data point. High sensitivity, 
high resolution data in the ran ge ^ 100 — 1000 G Hz, such 
as those to come from Planck jTauber et al.ll2010ll . are re- 
quired to rule this out and to obtain precise values for the 
thermal dust contribution. 

4.3 Anomalous microwave emission 

In Sections l4.ll andl 4.2l we used multi-frequency data to place 
limits on free-free and thermal dust emission in LDN1621 at 
31 GHz, as observed by the CBI. The bulk of the emis- 
sion at 31 GHz appears to be anomalous in that it can- 
not be easily explained by the usual diffuse components of 
the ISM at radio/microwave frequencies. Extrapolating the 
best-fitting thermal dust model curve to 31 GHz (Fig. [3]), 
gives a predicted 31 GHz flux density of 0.151 ± 0.079 Jy. 
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Table 2. Integrated flux densities for a 30 arcmin diameter aper- 
ture centred on LDN1621, at 13 arcmin resolution. *Note that 
the 31 GHz value was estimated by simulating the 31 GHz map 
based on the cross-correlation with the 100 /^m template. 



Frequency 


Flux density 


(GHz) 


(Jy) 


31* 


1.67 ±0.66 


93.5 


1.69 ±0.62 


2997 


368 ± 62 


4995 


155 ± 16 


11988 


47.8 ±3.3 


24975 


26.5 ±2.3 



1000.00 



10.00 r 




100 1000 
frequency (GHz) 



10000 



Figure 3. Spectrum of LDN1621 for 13 arcmin resolution data, 
within a 30 arcmin aperture (see text). The solid line represents 
the best-fitting modified blackbody curve using the 93.5, 2997 
and 4995 GHz data points. The predicted 31 GHz flux density is 
shown as an open square. The actual observed value (estimated 
via cross-correlation with the 100 /^m map; see text) is shown as 
a filled circle. The dashed line is for a fit where the emissivity 
index is fixed at /3 = +1.7. 



This is well below the estimated true value at 31 GHz of 
1.67 ±0.66 Jy. There is therefore an excess of 1.52 ±0.66 Jy, 
which is significant at the 2.3(t level. For a fixed emissiv- 
ity index, /3 = ±1.5, the significance reduces to l.Scr. For 
an extreme value of /3 = ±1.0, the flattest value that is ex- 
pected on theoretical grounds, the significance is just O.Scr. 
High frequency data [u > 100 GHz), such as those expected 
from the Planck satellite, will be crucial in determining the 
contribution from thermal dust emission. 

A number of possibilities exist to explain the ex- 
cess 31 GHz emission inclu ding electro-dipole radiati on 

from spinning du st gr ains (iDraine fc Lazari"anl Il998al lbl: 

lAli-Haimoud et al.ll2009l: Ysard et al.ll20ld ). magneto-dipole 
radiation ( Draine &: Lazarian 19991 ). free- free from hot T ~ 
10^ K gas (iLeitch et al 19971) . flat-spectrum (hard) syn- 



chrotron (|Bennett et al 



. . , ■ I2OO3I) and low-level solid-state 

structural transitions ( Joned |2009| ). With the current ra- 
dio/IR data, it is somewhat difficult to confidently distin- 
guish between these possibilities. In other clouds, and in the 
diffuse ISM, the most favoured explanation is in terms of 
spinning dust grains. 

Several authors have used correlation analyses 
against dust templates to try to isolate the type of 
dust grains that correlate best with the radio data 
(e.g. Ida Oliveira-Costa et all l2002l : ICasassus et"al] I2OO6I : 



I Ami Consortium et aLll2009l : lYsard et al.ll2010l ). It is clear 
from Fig. [2] that there is a good correlation with the IRIS 
maps at 12 — 100 /.im, suggesting a possible dust origin. 
We calculated Pearson's correlation coefficient, r, between 
the CBI 31 GHz map and IRIS data at 12, 25, 60 and 
100 fim. We take into account the spatial filtering of the 
CBI instrument by simulating CBI observations, based on 
the real u, v data with IRIS maps as inputs, and making 
CLEANed maps in the same way as for the CBI data. 
Errors were estimated by using Fisher's r to z transfor- 
mation and c alculating the 68 per cent confidence interval 
iFisher I (|l915l ). Since there is some concern regarding the 
primary beam correction beyond FWHM/2, we computed 
the correlation for different sized regions. Table [3] gives the 
r correlation coefficients for four regions. Given the noise 
level, and uncertainty in the primary beam correction in 
the outer regions, the 31 GHz-IR correlation is remarkably 
high with r ~ 0.6 — 0.8; significant at ^ lOcr. For spinning 
dust, we would expect a better correlation with shorter IR 
wavelengths (~ 10 — 30 /im) because these are dominated 
by the smaller (and generally warmer) dust grains which 
are responsible for the bulk of spinning dust emission 
lAli-Haimoud et al.|[2009l : lYsard et al.ll2010l ). Longer wave- 
lengths (> 100 /xm) are dominated by larger (and typically 
cooler) grains. We find no strong preference for a particular 
IR template or short versus long IR wavelengths. The 
exception is for the 10 arcmin region centred in the middle 
of the ring, where there is a slightly better correlation with 
short (12/25 /im) wavelengths (r = 0.82 ± 0.06) compared 
with longer (60/100 /xm) wavelengths (r = 0.69 ± 0.09). 
This is consistent with the expectation for spinning dust. 

A number of authors have calculated the coupling 
coefficient between IR templates and the radio data. 
The coefficient for 100 ^im, sometimes referred to as the 
emissivity, is ~ 10 /tK (MJy/sr)~^ in the diffuse ISM 
with a scatter of a factor of ~ 2.5 (|Bandav et al.l l2003l : 
iDavies et al] I2OO6I ). In LDN1622 the coupling coefficient 
is 21.3 ± 0.6 /iK (MJy/sr)~^ while HII regions are typi- 
cally at < 5 liK (MJy/s r)~^ (|Dickinson et al.l |2006| . [2Ci07l . 
l2009al ; [Scaife et aUliooa ). In LDN1621, taking no account 
of any free-free contribution, the coupling coefficient is 
16—22 /iK (MJy/sr)~^ depending on which part of the image 
is used. For a 10 arcmin radius circle centred in the middle 
of the ring (see Fig.[l]), where the spatial correlation is high- 
est (Table El, the coeflicient is 18.1 ± 4.0 /iK (MJy/sr)~\ 
Assuming there is no free-free contribution, which would re- 
duce the coefficient, LDN1621 appears to emit at the same 
level as LDN1622 compared to the 100 /im template, and 
a factor of ~ 2 more than the average value found in the 
diffuse ISM. 



5 CONCLUSIONS 

LDN1621 is a region of diffuse emission ^ 25 arcmin to the 
north of LDN1622. Observations with the CBI at 31 GHz 
show a broken ring of emission, that is strongly correlated 
with FIR emission at 12 — 100 /im, with Pearson correlation 
coefficients in the range f» 0.6 — 0.8. Optical and Ha data 
show absorption of a strong background of emission from 
warm ionized gas in the Eastern arm of Orion. This suggests 
that LDN1621 and LDN1622 are in the foreground of Orion 
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Table 3. Pearson correlation coefficients between the 31 GHz 
data and IRIS maps at 12, 25, 60 and 100 lira. These were cal- 
culated over circular regions with radii given in the first column. 
The regions are centred on the middle of the image except for 
the aperture marked with a asterisk (*) which was centred in the 
middle of the ring feature at R.A. 05^54™46=, Dec. +02° 14^03'' 
(see Fig. [TJ. 



Radius 


12 lira 


25 ura 


60 /^m 


100 fira 


(arcmin) 










10* 


o.78±i).oa 


0.85±"-»'' 






10 


0.75±0" 


0.70±oi2 


0.75±0-" 


0.7i±o.i2 


15 


0.58±oio 


0.63±0"9 




0.57±o.io 


20 


0.63±007 


0.68±0-0^ 







S at a distance of ~ 500 parsec), possibly as close as 120 pc 
Wilson et alll2005l ). No Hq emission, associated directly 
with LDN1621, is seen. This suggests that LDN1621 itself 
is not emitting significant free-free emission, although the 
effects of dust extinction do not allow a strong constraint 
to be placed. Low frequency radio data also do not show 
evidence of diffuse emission associated with LDN1621. The 
31 GHz emission is at « 20 — 30 mjy beam~^ while an 
analysis of the GB6 map at 4.85 GHz provides a strong 
(3cr) upper limit of 7.2 mJy beam~^ at 31 GHz for free-free 
emission. The FIR-correlated emission at 31 GHz therefore 
appears to be mostly due to radiation associated with dust. 

IRAS data alone do not allow a reliable extrapolation 
of the Rayleigh- Jeans thermal dust tail to 31 GHz. WMAP 
data at 93.5 GHz combined with IRAS data allowed the flux 
density to be estimated in an aperture of diameter 30 arcmin 
at an angular resolution of 13 arcmin. A single modified 
blackbody indicates that the thermal dust is ~ 10 per cent 
of the total 31 GHz flux, corresponding to an excess of 1.52± 
0.66 Jy (2.3(t). The dust-correlated emission has a coupling 
coefficient, relative to 100 ^ra, of 18.1±4.4 (MJy/sr)"^, 
consistent with that observed from LDN1622. 

Orion East (consisting of both LDN1621 and LDN1622) 
appear to be part of the same system of dust clouds, emitting 
significant anomalous emission at frequencies ~ 30 GHz. 
Spinning dust is an obvious candidate for the physical 
mechanism responsible for the bulk of the emission. High 
sensitivity data, covering a wide range of frequencies (~ 
5 — 300 GHz), is required to study such clouds in more de- 
tail. Data from the Planck satellite will be particularly useful 
in constraining the Rayleigh- Jeans dust tail, which may be 
responsible for a significant fraction of the 31 GHz if the 
emissivity index flattens at longer wavelengths. 
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